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1 Introduction 

Consider the equations of motion of Bohmian mechanics for a system of N particles 
with masses mi, ...,mjy moving in physical space H 3 : the wave function ip evolves 



according to Schrodinger's equation 



(-E^A fe + ^))^), (1) 



dt 

and the configuration Q = (Qi, . . . , Qjy) G ]R 37V , with G 1R 3 denoting the 
position of the k-th particle, evolves according to Bohm's equation 

^ = **W, (2) 
where the velocity field = (vf , . . . , v^) is determined by the wave function ?/> 

v fc (<?) = — Im — (3) 
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(Bohm 1952). This velocity field is regular on the subset of H 3 where ip 7^ and is 
differentiable. The question arises as to what happens if the configuration Q t , when 
moving in accordance with @ along integral curves of v** , reaches at some time r a 
singularity of v ^, for example a node of the wave function, i.e. a point where ip = 0. 
In general, the event of reaching a singularity of corresponds to a singularity of 
the motion: the velocity dQt/dt becoming infinite, Q t being discontinuous, or even 
"exploding," i.e. reaching infinity as t — > r. In some cases it may be possible to 
continue the Bohmian trajectory through a singularity of the velocity field, but in 
any case, Eqn. (Q) is not defined at singularities of xr , and the theory would have to 
be supplemented by a suitable prescription for how to extend the motion through a 
singularity, or which trajectories to avoid. 

We should therefore consider the following problem — the initial value problem in 
Bohmian mechanics: For given initial values ipo and Qo at a time to (we shall put 
to = 0), do there exist global unique solutions (ip t ,Qt) of such that ipf = ipo 
and Q to = Qo ? A positive answer to this question for all or at least for a suitable 
majority of initial conditions, or for a class of "physically relevant" initial conditions, 
is certainly important: the limits of global existence of solutions may hint at limits 
of validity of the theory. 

Let us first give a simple explicit example where indeed some trajectories reach 
nodes of the wave function. Consider the one-dimensional harmonic oscillator (with 
h — m — 10 — 1), i.e, the potential V(q) = q 2 /2, and take as the wave function a 
superposition of the ground state and the second excited state: 

= e~fl\\ + (1 - 2q 2 )e- 2U )e- u / 2 . (4) 



2 



This wave function leads to a velocity field which is an odd function of q and 
periodic in time with period tt. The Bohmian motion is therefore invariant under 
reflection q — > — q and periodic in t with period tt: If Q(t) is a solution of @, then 
—Q(t) is also a solution of (H), and Q(£ + 7r) = The wave function if> has nodes 

at (g, t) = (0, (n + l/2)7r) and at (q, t) = (±l,n7r) for all integers n. There are three 
trajectories which periodically run into nodes of if), see Figure [|. One of them, the 
constant Q t = for t ^ (n+ 1/2)tt, is a solution of (0) which runs with velocity at 
times t — > (n + l/2)7r into nodes of ^. The other two "node-crossing" trajectories 
(which are reflection images of each other) are singular at the nodes: for example, 
in the vicinity of the node at q — 1, t = the trajectory running into it has the 
form Q t = (3t 2 /4) 1 / 3 + 1, i.e. it has infinite velocity at t = 0. 

In this example the trajectories may be continued through the nodes in an ob- 
vious and consistent way. This is an artifact of the low dimensionality: In fact, in 
d — 1 dimension, Q t satisfies 

/ M<l)\ 2 dq= \Md)\ 2 dq, (5) 

employing the equivariance of the |^t| 2 -measure, i.e. that the density po = iV'ol 2 
evolves under the Bohmian dynamics to p t = \ip t \ 2 (Diirr, Goldstein, Zanghi 1992a; 
see also Section 4), and the property that trajectories do not cross in configuration- 
space-time. Clearly, much less regularity of if) is needed for this definition of the 
motion of the Bohmian configuration — in particular, nodes of the wave function are 
no problem — and, if both definitions (0) and (|) are possible, they agree. However, 
a generalization of (|5|) to the physically most interesting case of d = 3iV-dimensional 
configuration space is not known. We shall turn now to the question of the existence 
of global unique solutions of the equations of motion of Bohmian mechanics, (P and 
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2 The Schrodinger equation 

In Bohmian mechanics, the evolution of the wave function ip according to 
Schrodinger 's equation (P is independent of the evolution of the actual particle 
configuration Q according to the guiding equation (|2]), while for the integration of 
(|2p we need ip. Therefore, when solving the initial value problem, we may consider 
Schrodinger's equation first. 

The linear partial differential equation (|I|) is usually discussed in the frame- 
work of the Hilbert space Tt = L 2 (H 3Ar ) of square integrable functions, by viewing 
~ N h 2 

the Hamiltonian H = — ^ + V(q) as the generator of the unitary group 

k—l k 

U t = e~ %m l n (Kato 1951, Reed and Simon 1975). Usually, boundary conditions 
on the wave function have to be specified in order to get a unique time evolution 
{Ut)tem- m more technical terms, one has to select a self-adjoint extension H of the 
partial differential operator H, which is a priori defined only on sufficiently smooth 
functions. As different self-adjoint extensions generate different time evolutions of 
the wave function, the choice of the right self-adjoint extension is a matter of physics. 
We shall comment on this in Section 4. Given (Ut)teiRi f° r an y initial ipo £ "H we have 
a global unique wave function ip t = Utipo- This wave function, however, in general 
is not a genuine solution of Schrodinger's equation: Generic ip EH are not smooth 
functions, and ipt will not be differentiable, so that the differential equation (|J) 
cannot be discussed. The standard procedure now is to forget about Schrodinger's 
equation, and argue that all one is interested in is a unitary time evolution — or, 
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even more abstractly, a unitary representation of the time translation group. This 
is, however, not sufficient for our approach. From the point of view of Bohmian 
mechanics, the wave function ip is a smooth field on configuration-space-time solv- 
ing Schrodinger's equation. To obtain this, we have to put suitable conditions on 
the initial wave function ip . It turns out that a suitable subspace of H is given 

oo 

by the so-called C^-vectors of the Hamiltonian H, C°°(H) = f) V(H n ). V(H n ) 

n=i 

denotes the domain of the n-th power of the Hamiltonian, i.e. the set of wave func- 
tions for which the expectation value of the 2n-th power of H, the "energy," is finite, 
/ 3N (H n ^) 2 dq = (^\H 2n ^) < oo. 1 C°°(H) is a dense subset of H, and invariant un- 
der the time evolution U t . Eigenfunctions and wave packets are special C 00 - vectors, 
so all wave functions usually considered in physics are included. It should not at all 
be regarded as a defect of Bohmian mechanics that it is not defined for all tpo £ 'H- 
From the point of view of Bohmian mechanics, the Hilbert space L 2 (JR 3N ) is not 
the state space of the wave function, but a useful tool for the analysis of the theory. 
(In this context, see also the contributions of Diirr, Goldstein, Zanghi and Daumer 
to this volume.) 

We have: For any ij e C°°(H), ip t = U t ip = e^/Vo is a global smooth solu- 
tion of Schrodinger's equation on Qx IR, where Q denotes the set where the potential 
is smooth. This result cannot be genuinely strengthened: the wave function cannot 
be expected to be regular at points where the potential is singular. For instance, 
the ground state eigenfunction of the Coulomb potential V(q) = l/\q\, e - ' 9 ', is not 
differentiable at the singularity of the potential at q — 0. Important examples of 
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potentials are the iV-particle Coulomb interaction for particles with charges 

N N CiCj 

Vbouiomb(qi) • • • , qjv) = X! X! I — zr — r 

i=l j=i+l 11* Oil 

with ^coulomb = IR 3Ar \ U U "tl* = %}> an< ^ ^ ne -^- e l ec t ron atom 

i=l j=i+l 

N N N 

atom(qi, . . . , qjv) = 2^ TT + 2-2- I — Zi 

i=i I9.il i=ij=i+i lli Oil 

/AT N N \ 

with fi a tom = lR 3Af \ Ui^* = 0}U (J [J {qj = , in the approximation that 

Vi=l i=lj'=t+l / 

the nucleus with charge e n is at rest at the origin, acting like an external Coulomb 
field. 



3 The Bohmian trajectories 

With the wave function ip G C°°(Q x H), the velocity field ([3D can be formed on 
the set where the wave function ip ^ 0. We shall call the set where the velocity field 
is regular Q = (Q x 1R) \ A/", where M is the (space-time) set of nodes of the wave 
function M = {(q,t) G Q x 1R : ipt{q) — 0}. On this set, the first order ordinary 
differential equation (0) is locally integrable. If we extend the solution as far as 
possible, we obtain for each initial value go £ Go, i-e. (go, 0) G a maximal solution 
Q(t;qo) on a maximal time interval of existence {t~ (qo) , t + (q Q )) . The solution is 
called global if t~ = — oo and r + = +oo. In the introduction we have given an 
example showing that we cannot expect to have global solutions for all initial values 
go- But what indeed holds true is that for almost all initial configurations we 
have global unique solutions. "Almost all" here is with respect to the natural — the 
equivariant — measure on configuration space H 37V , that is the measure with the 
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density (V'ol 2 - We have proved that for a large class of potentials, including the N- 
particle Coulomb interaction with arbitrary masses and charges, as well as arbitrary 
positive potentials, P^°(r + = +00 and t~ = —00) = 1. In other words, for almost 
any initial point Q , the solution of the guiding equation (§) is global, i.e. nodes 
of the wave function or other singularities of the velocity field will not be reached 
in finite time, and the solution does not "explode," i.e. reach infinity in finite time. 
(Cf. the example of Figure [1|: the set of initial configurations for which the solution 
of (0) is not global consists of 3 points — certainly a set of P^°-measure — while for 
all the other initial values the solution is global.) 

We shall comment on the proof — which is in fact quite intuitive — in the next sec- 
tion. (For the details, see Berndl, Diirr, Goldstein, Peruzzi, Zanghi 1995.) Despite 
the relative simplicity of the proof, the generality of the result is rather surprising. 
The analogous problem in Newtonian mechanics — global existence and uniqueness 
of solutions of the iV-body problem for Newtonian gravity — is a classical problem of 
mathematical physics which has been investigated with a great variety of methods 
(Moser 1973, Diacu 1992), but for which the analogous result — global existence of 
solutions for Lebesgue-almost all initial values in the iV-particle phase space H 6Ar — 
has not yet been established. In addition to the possibility of collision singularities, 
the iV-body problem with N > 3 yields marvelous scenarios of so-called pseudocol- 
lisions, where some particles, while oscillating wildly, reach infinity in finite time. 
Examples of such catastrophies have been constructed by Mather and McGehee 
1975, 2 by Gerver 1991, and by Xia 1992. For special situations, such as for example 
an almost planar solar system with weakly eccentric planets, the KAM theorem 
furnishes among other things the existence of global solutions for a set of initial 
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values of positive Lebesgue measure in phase space (cf. for example Arnold 1963). 
For the general case of the iV-body problem with N > 5 it is still open whether for 
almost all initial values the solutions of Newtonian gravity are global or whether 
pseudocollisions occur on a set of positive measure. 

The explanatory power of Newtonian mechanics rests largely on the analysis 
of specific solutions of the equations of motion, such as for example the planetary 
motion, or the gyroscope motion — physical situations where the bodies are modelled 
by fixed mass densities. The more ambitious program of Newtonian mechanics for 
a system composed of a huge number of point particles or rigid balls, for example 
the statistical gas theory, is difficult for two reasons. Firstly, to get global existence 
of solutions for almost all initial values — certainly a prerequisite of a statistical 
analysis — one has to deal for example with the problem of singular interactions in 
the case of point particles resp. multiple collisions in the case of rigid balls. Secondly, 
more interesting than equilibrium properties and much more difficult to analyse are 
the physical effects occuring in the transition from nonequilibrium to equilibrium. 

In Bohmian mechanics, actual trajectories are interesting in some cases. How- 
ever, to determine the equilibrium properties they are not needed. And in contrast 
to the situation in classical physics, we have the strongest empirical evidence that 
our world is in quantum equilibrium (Diirr, Goldstein, Zanghi 1992a and 1992b). 
Our theorem of global existence of Bohmian trajectories is thus exactly what is 
necessary and sufficient for a complete analysis of Bohmian mechanics. It provides 
the rigorous basis for the derivation of the quantum formalism as well as scattering 
theory from an equilibrium analysis of Bohmian mechanics. (These topics are de- 
scribed in the contributions of Diirr, Goldstein, Zanghi and Daumer to this volume.) 
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In this sense, what we have proved is just that \ip\ 2 is an equivariant density for the 
Bohmian particle motion, 

Po = |^o| 2 Pt=\A? foralUeR. (6) 

4 The quantum flux 

The preceding sentence may have confused readers: Why isn't the equivariance of 

the |^| 2 -measure clear without an intricate examination of the global existence and 

uniqueness of solutions of Bohmian mechanics? Why doesn't © follow immediately 

from a comparison of the continuity equation for an ensemble of configurations 

moving with velocity and having a density pt{q) 
r) N 

%-pt(q) + £div fc (vf(g)p t (g)) = Q (7) 
01 k=i 

with the "quantum continuity equation" 

a N 

q:\Mq)\ 2 + Ediv,jf(?) = 0, (8) 

k=l 

noting that the quantum probability current = (jf , ■■-,]%) is given by 

f k = vfl^l 2 = -Im(fV^) ? 
m k 

The answer is that the continuity equation (|7|) expresses how an ensemble density 
Po evolves under the deterministic evolution of the trajectories, and holds therefore 
on just that set which is covered by the integral curves of v^. We shall denote this 
set — the image set of the maximal solution Q — by I: 

I = {(q,t) E Q : 3q G Go with t E (r _ (q ) , r + (q )) and q = Q t (qo)}- 

Equation (|3|), on the other hand, is an identity for every ip t which solves Schrodinger's 
equation. Thus © holds on I: pt(q) = |'0t(?)| 2 f° r fat) G I- By putting ptfa) = 
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on the set which is not reached by the Bohmian trajectories (Q x H) \ X, we obtain 
Pt(q) < Itftio) | 2 f° r (?> t) £ f2 X M. This insight is fundamental to the proof of global 
existence. 

The second fundamental insight comes from the consideration of the 
configuration-space-time flux J(q,t) = (v^ t (q)pt(q), Pt(q))- To establish P^o(Y+ = 
+00 and r~ = —00) = 1, i.e. that for P^°-almost all initial configurations go the 
maximal solution is global, we have to show that the solution will reach singularities 
of the velocity field or infinity in finite time at most for a set of initial configurations 
of P^°-measure zero. The probability of such a "bad event" is estimated with the 
help of the flux J(q,t): For any hypersurface £ in configuration-space-time with 
local normal vector field n(q, t) and surface area element da, 



JjJ(q,t)-n(q,t)\da 



is the expected number of crossings of the hypersurface £ by the Bohmian trajectory 
and hence a bound for the probability of crossing E. We obtain 

(Q t crosses S) < J \ J(q, t) ■ n(q, t)\da< J \ J* (q) ■ n(q } t) | da (9) 

with the quantum flux J^(q) = (j^(q), \ipt(q)\ 2 )- 

Consider now neighborhoods around the singularities of the velocity field: A/" e , a 
(configuration-space-time) neighborhood of thickness e around the set of nodes M of 
the wave function, S s , a (configuration space) neighborhood of thickness 6 around 
the set of singularities of the potential S = dQ, and JC r , a sphere in configuration 
space of radius r to control escape to infinity. Q e5r denotes the set of "e-5-r-good" 
points in configuration-space-time: Q eSr = ((/C r \ S s ) x H) \ Af e (see Figure 0). 
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A Bohmian trajectory approaching a singularity of the velocity field or infinity 
first has to cross the boundary of Q eSr . From @, we obtain the following bound for 
the probability of a "bad event" in the time interval [0, T\. for all < T, e, 5, r < oo 

P*»(t+<T) < P^(G \Gt) + I t 3N \J^(q)-n(q,t)\da 

Jag €5r n(u 3N x[o,t]) 

< P^(Go\Go r ) + N(e,«y,r) + S(5) + I(r) (10) 



with 



N(e,5,r) = [ \J*(q) ■ n(q,t)\da, 

S(6) = I \J^(q)-n(q,t)\da, 

I(r) = / \J^(q) -n(q,t)\da. 

J(dlC r nU)x\o,T} 



'(9x; r nn)x[o,r] 

By proving that for appropriate choices of sequences e — >• 0, 5 — > 0, and r — > oo the 
right hand side of ( jlOl) gets arbitrarily small we show that for all T > 0, P^°(r + < 
T) = 0, and thus P^°{r + < oo) = 0. From time reversal invariance we obtain that 
also P^°(r _ > — oo) = 0, and that altogether the solutions of Bohm's equation (0) 
with the velocity field ((^) are global for almost all initial configurations. 

Heuristically, it is rather immediate that the flux integrals N, S, and I get 
arbitrarily small as e — > 0, 5 —>■ 0, and r — > oo. For N, observe that the flux = 
at the nodes of the wave function. From the continuity of ip, is small in the 
vicinity of Af. Furthermore, the nodal set TV itself is generically small: since if} is 
a complex function, the set where ipt{q) = 0, i.e. Keif)t(q) = and lmif>t(q) = 0, 
is generically of codimension 2 in configuration-space-time. Thus the area of dM e 
should be small. Of the set S of singularities of the potential V, we assume in 
Berndl, Diirr, Goldstein, Peruzzi, Zanghi 1995 that it be contained in a union of a 
finite number of (3iV — 3)-dimensional hyperplanes, as is certainly the case for the 
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iV-particle Coulomb interaction V = Vbouiomb and the iV-electron atom V = Ktom- 
The area of dS s is therefore very small. Moreover, in certain cases it is required as 
a boundary condition for self-adjointness of the Hamiltonian that the flux into the 
singularities of V vanishes, see below. That the flux to infinity is small can be derived 



follows from the square integrability of ip and V?/>, i.e. from the normalizability of 
the |^>| 2 -distribution and finite "kinetic energy." These conditions are automatically 
fulfilled for the considered class of potentials for tp e C°°(-f/'), which is just what 
we required for the existence of a global classical solution of Schrodinger's equation. 

We are thus lead back to the question of the existence of global classical solu- 
tions of Schrodinger's equation! In fact, the connection between the global existence 
of Bohmian trajectories and global solutions of Schrodinger's equation is most re- 
markable. The clue is the quantum flux J^, which in Bohmian mechanics has the 
interpretation of a flux of particles moving along deterministic trajectories with ve- 
locity The condition that there is no flux into the critical points ensures firstly, 
as explained above, that the Bohmian configuration will not reach the critical points 
and thus exists globally. Secondly, it provides suitable boundary conditions for the 
domain of the Hamiltonian T>(H) such that the Hamiltonian will be self-adjoint on 
T>{H) and thus Schrodinger's equation has global unique solutions as explained in 
Section 2. This connection is realized by considering that from integration by parts 



from the fact that the quantum flux J^*{q) tends rapidly to as \q\ — > oo, which 




for M = IC r \S s , and thus by the self-adjointness of the Hamiltonian 
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This is only slightly weaker than the vanishing of S(S) and I(r) in the limit 5^0, 
r — > oo, which is part of our sufficient condition for global existence of Bohmian tra- 
jectories. Moreover, in situations where the self-adjoint extension of H (cf. Section 
2) is not unique, the particle picture of Bohmian mechanics supplies an interpreta- 
tion of the different possible boundary conditions yielding different time evolutions 
of the wave function and thus a basis for the choice of one over the others. (For 
more on these points, see Berndl, Diirr, Goldstein, Peruzzi, Zanghi 1995.) 

In this way, the point of view of Bohmian mechanics provides genuine under- 
standing to the mathematics around the self-adjointness of Schrodinger Hamilto- 
nians! Usually, the self-adjointness of the Hamiltonian — via its equivalence to the 
existence of a unitary group — is motivated by the conservation of |^| 2 -probability. 
Probability of what? The standard answer — the probability of finding a particle in 
a certain region — is justified by Bohmian mechanics: A particle is found in a certain 
region because, in fact, it's there. By incorporating the positions of the particles 
into the theory, and thus by interpreting \ip\ 2 as a probability density of particles 
being and the quantum flux as a flux of particles moving, Bohmian mechanics 
can be regarded as providing the foundation for all intuitive reasoning in quantum 
mechanics. 
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Figure 1: Sketch of the course of Bohmian trajectories for the wave function (4). 
Nodes of the wave function are marked by dots. 
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Figure 2: Bohmian trajectories run in the white area Q e&r = ((/C r \ S 5 ) x IR) \ Af e . 
The set of singularities of the potential S x IR is dashed, nodes of the wave function 
are marked by x. Trajectories having crossed dQ e5r are dotted. By letting the grey 
neighborhoods Af e and S 5 around the singularities of the velocity field shrink and 
the sphere K7 blow up, the set of dotted trajectories shrinks to a set of measure 
zero. 
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Footnotes 



1 If only the expectation of H 2 , the squared energy, is finite, vp G V(H), 
"Schrodinger's equation holds in the L 2 -sense," i.e. ^lim — — = Htp t , where 
the equality and convergence of the limit holds with respect to the Hilbert space 
norm and not pointwise as expressed by Schrodinger's equation. 

2 In this (one-dimensional) example the system explodes only after infinitely 
many two particle collisions; thus it does not describe a genuine pseudocollision. 
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